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We present a new, simple method for calculating the scalar, electromagnetic, and gravitational self 
forces acting on particles in orbit around a Kerr black hole. The standard "mode-sum regularization" 
approach for self-force calculations relies on a decomposition of the full (retarded) perturbation 
field into multipole modes, followed by the application of a certain mode-by-mode regularization 
procedure. In recent years several groups have developed numerical codes for calculating black 
hole perturbations directly in 2+1 dimensions (i.e., decomposing the azimuthal dependence into m- 
modes, but refraining from a full multipole decomposition). Here we formulate a practical scheme 
for constructing the self force directly from the 2+1-dimensional m- modes. While the standard 
| mode-sum method is serving well in calculations of the self force in Schwarzschild geometry, the 

, new scheme should allow a more efficient treatment of the Kerr problem. 
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I. INTRODUCTION 



The motion of a small particle (a point mass, a point electric charge, or a point scalar charge) in curved spacetime, 
in situations where the perturbation caused by the particle can be considered "small" , may be described in terms of an 
effective self force (SF). The SF accelerates the particle with respect to the "background" spacetime; it contains both 
conservative and dissipative parts, the latter being interpreted as the radiation reaction force. A formal expression 
for the electromagnetic (EM) SF in curved spacetime was derived in 1960 by DeWitt and Brehme [l| (later corrected 
by Hobbs [2]). In 1997, Mino, Sasaki and Tanaka Q derived a formal expression for the gravitational SF, in two 
different ways: by extending the method of DeWitt and Brchcmc, and also using the technique of matched asymptotic 
expansions. An alternative derivation of both EM and gravitational SFs was presented by Quinn and Wald [J]. The 
formal expression for the scalar SF was obtained by Quinn [5[ in 2000. Detweiler and Whiting |6j later showed how the 
SF interpretation of the motion is consistent with the standard picture of geodesic motion in a perturbed geometry. 
A thorough self-contained review of these developments can be found in Q • 

For all fields (scalar, EM, or gravitational), the formal expression for the SF contains a "tail" term, which describes 
the interaction of the particle with the part of its field supported inside the particle's past light-cone. In the grav- 
itational case, the tail part of the self-interaction is responsible for the entire SF effect (for particles moving along 
1 momentarily-geodesic trajectories in vacuum spacetimes). In practice, to facilitate actual calculations of the SF in 
\ concrete situations, it is convenient to express the tail force as the difference between the "full" force, arising from the 
full (retarded) perturbation field, and the "direct" force, describing the back reaction from the instantaneous part of 
| the retarded field supported only on the particle's past light-cone. 

The implementation of the above theoretical framework in the concrete problem of a particle in orbit around a 
black hole has been the subject of intensive study over the last decade. (This has been largely motivated by the need 
to accurately model the orbital evolution of astrophysical compact objects inspiralling into massive black holes — 
of the prime targets for LISA, the planned space-based gravitational wave detector.) To calculate the SF in this 
scenario, one usually starts by solving (numerically) the appropriate linear perturbation equations over the black hole 
background, with a stress-energy source corresponding to the orbiting particle. One normally exploits the symmetry 
of the black hole background, by writing the perturbation equations in a separated form and tackling them mode by 
mode. The main practical challenge, then, is to correctly split each of full, retarded-field modes into its "tail" and 
"direct" contributions. A scheme to achieve this was first devised by Barack and Ori in 2000 [8|, and later developed 
to deal with arbitrary geodesic orbits in Schwarzschild d E3, EH EH > anci , eventually, arbitrary geodesic orbits in 
Kerr [l3j]. This "mode-sum regularization" scheme has since been implemented by various authors on a case-by case 
basis, so far only for orbits in Schwarzschild [3, EM EM [H El, El Hfl HI- Most recently, the mode -sum scheme was 
used to calculate the scalar SF for generic geodesic orbits in Schwarzschild [12] , and the gravitational SF for circular 
geodesies in Schwarzschild [23j j . 

The standard mode-sum scheme uses as input the individual multipole modes of the full perturbation field (which 
are normally obtained numerically). The contribution of each mode to the full force (which is finite at the location of 
the particle) is subjected to a certain regularization procedure, which essentially amounts to subtracting out the direct 
force mode-by-mode. The sum of the regularized modes is guaranteed to converge and yield the correct tail force. 
In this description, the term "multipole mode" refers to a particular spherical-harmonic component I of the field, 
summed over all azimuthal numbers m. In spherically-symmetric black hole spacetimes (including Schwarzschild), 
these £-modes can be obtained using either frequency-domain or time-domain methods. In the former approach (e.g., 
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[20|) one first solves the ODE for each Fourier frequency to of each multipole £, m of the field, and then sums the 
contributions from all u> and m for given I. In the second, time-domain approach (e.g., [H,[23j]), one solves a PDE in 
1+1D (radius+time) for each £, m, and then sums over m. 

The situation in Kerr is slightly more involved, as in this case the field equations do not separate into individual £, m 
modes in the time domain, and the 1+1D approach is not applicable. In this case, the ^-modes required as input for 
the mode-sum scheme are obtained by first solving the ODE for each (spin-weighted) spheroidal-harmonic £, to, lu of 
the perturbation, then summing over u>, to for given i, and finally re-decomposing each spheroidal-harmonic ^-mode in 
a basis of spherical harmonics I. Although this procedure is mathematically quite straightforward, its implementation 
may be rather cumbersome. It has not been attempted so far. 

An alternative, more direct approach in the Kerr case may become possible thanks to recent developments which 
facilitate the calculation of black hole perturbations in 2+1D. In this calculations one exploits only the axial symmetry 
of the underlying Kerr geometry: The field equations are decomposed only into azimutal modes (each oc e mtp , where 
tp is the Boyer-Lindquist azimuthal coordinate), and one solves directly for the 2+1D TO-modes using time evolution. 
Over the past decade, several authors have considered evolution in 2+1D, with or without a particle source. Krivan et 
al. [lil, HH explored the late-time decay of scalar fields and Weyl-scalar vacuum perturbations by evolving the master 
Teukolsky equation in 2+1D. Pazos-Avalos and Lousto 2G] presented an improved, fourth-order-convergent code in 
2+1D, for the evolution of vacuum perturbations of the Teukolsky equation. Particle orbits in Kerr were tackled 
in 2+ ID by Lopez- Aleman et al. (27|, Khanna [28j and Burko and Khanna [29| . using a Gaussian representation 
of the particle. A more sophisticated finite-impulse representation was very recently suggested and implemeted by 
Sundararajan et al. [3(J. Sopuerta et al. proposed the use of finite-element methods for an effective treatment of the 
particle in 2+1D. This idea was implemented so far only in a 1+1D context [3l|,[32], and it awaits further development. 
Most recently, Barack and Golbourn [33| proposed a "puncture" scheme (further discussed below) for dealing with 
the singular behavior of the field in 2+1D. They demonstrated the applicability of this method in the test case of a 
scalar charge set in a circular orbit around a Schwarzschild black hole. 

Suppose that, for a given orbital configuration in Kerr, we had at hand the 2+1D m-modes of the retarded 
perturbation field (say, in the form of numerical solutions). We could then proceed by decomposing these modes into 
spherical-harmonic £, TO-modes, which, upon summation over to for given £ and evaluation of the contribution from 
each £-mode to the full force, would provide the necessary input for the standard mode-sum scheme. However, it would 
clearly be much more straightforward to construct the SF directly from the m-modes, avoiding the ^-decomposition 
altogether. The goal of this work is to formulate such an "m-mode regularization" scheme, which allows access to the 
SF directly from the 2+1D m-modes. 

There is an important difference between ^-mode and TO-mode formulations: The ^-mode of the full, retarded field 
(summed over to for given £) is continuous at the particle's location, and its contribution to the full force is finite. This 
is true for all types of perturbations, scalar, EM and gravitational (in the Lorenz gauge) alike. In contrast, the 2+ ID 
TO-mode perturbation diverges at the particle (logarithmically, in general [HI), and its contribution to the full force 
is divergent as well. This is troublesome in two ways: Firstly, it complicates significantly the numerical treatment. A 
numerical evolution scheme formulated in 2+1D, with a pointlike source, has to deal somehow with the divergence 
of the solutions at the particle's location. Secondly, since the individual m-mode contributions to the full force are 
divergent, they cannot be used, as they are, as input for an m-mode regularization scheme — at least not in the same 
way that the (finite, well defined) ^-modes are used as input for the standard mode-sum scheme. 

The puncture method of Ref. [33| offers solutions for both above difficulties. To illustrate the essential technique, 
consider the case of scalar perturbations (the EM and gravitational cases will be described later). Let <!> be the 
retarded perturbation field caused by a scalar charge in orbit around a Kerr black hole. In the puncture scheme we 
formally split <E> as 

$ = $ ros + $p, (1) 

where the "puncture field" $p is a certain function, given analytically, which approximates the singular behavior of $ 
well enough to guarantee that the m-modes of the "residual field" <& res are continuous at the particle's location [4l| . 
We then utilize $ res as a new variable for the numerical evolution: We re-write the field equation in terms of <!> r es (it 
will now have a source term which depends on the known function $p), separate the azimuthal dependence, and use 
evolution in 2+ ID to solve for each of the m-modes of $ rcs , which are continuous fields. The m-mode of the actual 
(divergent) field $ is then just the sum of the m-mode of $ res , obtained numerically, and the m-mode of $p, given 
analytically. This scheme deals with the first of the aforementioned difficulties. It does so, essentially, by treating the 
singularity of the field analytically, and solving numerically only for the residual, continuous field. 

The analysis in Ref. [33| focused on the calculation of the (scalar) field itself, and was not concerned with its 
derivatives. It incorporated a "leading-order" puncture, for which the m-modes of $ res were continuous but generally 
not differentiable (their derivatives diverged logaritmically at the particle). In order to deal with the second of the 
aforementioned complications — the divergence of the m-mode contributions to the full force — we will consider in the 
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present work an improved version of the puncture scheme, in which the m-modes of $ rcs are not only continuous, but 
also have continuous derivatives at the particle. We will prescribe a puncture function $p which achieves this. Our 
new mode-sum scheme for the SF will require as input the m- modes of the residual function $ res . We will show that 
the correct tail part of the scalar SF is simply the sum over all m-mode forces exerted on the particle by the m-modes 
of $ros, with no further regularization required. Similar results will apply in the EM and gravitational cases. 

This paper is structured as follows. Sections II, III, and IV are each devoted to one of the field types: Scalar, EM, 
and gravitational, in order. Each of these sections contains three parts. In the first part we review the formulation 
of the equation of motion with a SF term; in the second part we prescribe our m-mode regularization scheme for the 
SF; and in the third part we justify the scheme and explain the derivation of the new mode-sum formula. Section V 
contains a summary and some concluding remarks. 

Throughout this work we use standard geometrized units (with c = G = 1) and metric signature ( — | — h+). The 
Riemann tensor is defined as in Ref. [34| . and t, r, 9, ip are the standard Boyer-Lindquist coordinates. 



II. SCALAR FIELD 



A. Preliminaries 



Consider a test particle of a scalar charge q, moving freely in the vacuum exterior of a Kerr black hole with mass 
M ^S> q and angular momentum aM. Neglecting SF effects, the particle moves along a geodesic x^ = z^(t) of the 
background spacetime, parameterized by r, the proper time. We denote the full (retarded) scalar field associated 
with this particle by $(x), and assume that it satisfies the minimally-coupled Klein-Gordon equation, 

V Q V Q $ = -4vrp. (2) 

Here V Q denotes covariant differentiation with respect to the background Kerr geometry, and indices are raised and 
lowered using the background Kerr metric g a p. The scalar charge density on the right-hand side (RHS) is given by 

/>oo 

p{x) = q S^x-ziT^-gy^dT, (3) 



where g is the determinant of g a p, and x, z are short-hand for . 

Including self interaction of order q 2 (and ignoring the gravitational SF), the equation of motion of the particle can 
be written in the form [j| 

Vv^mh") = F« lf = lim F&y(x), (4) 

x — >z[t) 

where /i is the rest mass of the particle (assumed much smaller than M), u a = dz a /dr is the particle's four- velocity, 
and the "tail force" field is defined by 

F t Z il (x)=q 2 lim f \ a G rct [x,z( T ')]d T '. (5) 

Here G le t[x,x'] is the retarded Green's function associated with Eq. ([2|), V Q acts on the first argument of G re t, and 
T-(x) is the value of t at which the past light-cone of point x intersects the particle's worldline. Notice that in Eq. ^ 
we have kept /i inside the derivative. This is necessary, as \x will generally be time-dependent (this relates to the fact 
that, in general, the orbiting scalar particle will emit monopole radiation, a process which alters the particle's rest 
mass). The components of Eq. (0| orthogonal and tangent to u a give, respectively, the acceleration of the particle 
and the rate of change of ^.[42j| 

We comment on the regularity of the tail term, which will play an important role in our analysis: The limit in 
Eq. ([5]) cuts the worldline integral short of the light-cone singularity of the integrand at z = z(t_). As a result, the 
integrand is a smooth function of x even at x — > z. Since t^(x) is a continuous function of x [as x — > z(t) we have 
r_(x) — > t], the field F^^x) is continuous near x — z, and the tail term F' t " il [z(T)] in Eq. ^ is well defined. However, 
the integration limit t_(x) does not depend smoothly on the coordinates x, which impairs the smoothness of F^^x). 
More precisely, V/3T_(a;) has a finite-jump discontinuity at x ~ z [see Eq. (16) of [5J], which renders the derivatives 
of F^^x) discontinuous, yet bounded, at x = z. Our basic working assumption in the formulation of the m-mode 
scheme below will be that the tail-force field F^^x) is continuous for all x, and has at least piecewise continuous 
derivatives. 
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B. m-mode scheme: A prescription 



In the following we let z(r) be an arbitrary bound geodesic orbit of a scalar particle around a Kerr black hole, 
and prescribe the construction of the SF at an arbitrary point along this orbit, using the proposed m-mode scheme. 
The basic SF construction formula is Eq. (fT3"]) below. In this subsection we merely state this formula; in the next 
subsection we will explain its derivation. 

• Step 1: Construct the puncture function. 

For an arbitrary spacetime point x outside the black hole, let S be the spatial hypersurface i=const containing 
x, let f(t) be the value of r at which the particle's worldline is intersected by £, and denote z(t) = z[f(t)]. 
Define the coordinate distance 5x a = x a — z a (t), and construct the two quantities 

So = (g a p +u a u /3 )\ s Sx a Sx 13 , Si = (uAU 7 r^ + g a f 3n /2)\_8x a Sx l3 Sx 1 , (6) 

where the four-velocity u a , the Kerr metric g a /3, its derivatives g a p^ and the connection coefficients T^g are all 
evaluated at the worldline point z. Finally, define the puncture function as 



$ P = — , where e P = y/S + S%. (7) 
ep 

For Sx 7 very small we have S = Sq + Si + 0(6x ), where S is the squared geodesic distance from x to the 
particle's worldline (i.e., the squared length of the small spatial geodesic section connecting x to the worldline 
and normal to it). The above two leading terms of 5*, So(oc fa 2 ) and S\(oc Sx 3 ), are derived, e.g., in Appendix 
A of Ref. [Til] - We emphasize, however, that here the definitions in Eqs. ([6]) and (0 apply for arbitrary Sx, not 
necessarily small. 

• Step 2: Write down the field equation for $ rcs and separate into m-modes. 
Define the residual field 

$ rcs = $ - $ P (8) 

as in Eq. (fTJ), and re- write the scalar field equation {2J in the form 

V Q V Q $ rcs = ~4np - V Q V Q $ P = Z rcs . (9) 

The source Z rcs is extended (not confined to the particle's worldline), but contains no Dirac-delta on the 
worldline. As we show later, the field ^rcs is continuous at the particle, and its derivatives there are bounded 
(albeit generally discontinuous). 

Now formally decompose $ rcs and Z rcs into azimuthal m-modes, in the form 

oo oo 

$ rcs = Yl Cs(«, Z ICS = Ks(t,r,e)e^, (10) 



m— — oo m— — oo 



and use these expansions to separate the tp dependence in Eq. |9]). Each of the (complex- valued) m-modes 
4>™ s (t,r,6) will satisfy a hyperbolic field equation in 2+1D, of the form 

4 3) Cs = (ii) 

where Dg is a certain (m-dependent) second-order differential operator. The source modes are given explicitly 



by [13 



ZZ = ±f" Z^e-^'dtp', (12) 



which can be evaluated either analytically (as in [33]) or numerically. As we show later, the m-modes 
continuous and differentiable (have continuous first derivatives) at the particle. 
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• Step 3: Formulate an initial/boundary-condition problem for (f>™ s - 

We are looking for particular solutions (/>™ s which, through 0™ s + (j)™ = m , give the physical, retarded-held 
modes 0™. [Here 0™ and 0™ are the m-modes of $p and $, defined in analogy to (/>™ s in Eq. (p~0|) ] . For this 
to happen, the boundary conditions for </>™ s should be such that </>™ s + <fip represents purely outgoing radiation 
at the far "wave zone" (t, r 3> M), and purely ingoing radiation near the event horizon. Since the puncture 
modes </>p are given analytically (at least in terms of closed- form definite integrals), this readily translates to 
physical boundary conditions for 0™ s . In practice, it may be more convenient to re-define the puncture 4>p by 
(smoothly) suppressing its support away from the particle. The new residual function (j> les then coincides with 
(j) away from the particle, and the usual ingoing/outgoind boundary conditions will apply to it. A similar idea 
was implemented in the analysis of Ref . [331 ] . 

The choice of correct initial conditions for the numerical evolution is less of a concern. Linear perturbations 
on black hole backgrounds can be evolved stably over an indefinite amount of time. If the initial data for the 
evolution are specified much earlier than the time along the orbit where we wish to determine the SF, then the 
exact form of initial data (chosen reasonably enough) will have little effect on the outcome of our calculation. 
This, too, is demonstrated in the analysis of Ref. [331 ]. 

• Step 4- Solve for 0™ s using numerical evolution in 2+1D. 

We envisage solving Eq. (fTTj) using a finite-difference time-evolution code formulated in 2+1D, as in Ref. [33j |. 
Here we shall not be further concerned with the details of the numerical implementation, and proceed by 
assuming that the modes 0™ s have been calculated. 

• Step 5: Apply the m-mode-sum formula. 
The SF is simply given by 

oo 

w)i=«i;v a & , as) 

where (f>™ s (t,r,9,(p) are real fields constructed from cj)™ s (t,r,d) through 

C. = 2Re (C s e mv ) for m > 0, and Cf° = Cs = °- (14) 
As we argue below, the sum in Eq. (fT3"| is expected to converge at least as ~ 1/m. 

C. Analysis 

In what follows we justify and explain the above mode-sum prescription. We start by giving a detailed derivation 
of the mode-sum formula fj 13|) : we then explain its predicted convergence rate; and, finally, we analyze the behavior 
of the numerical integration variables 0™ s near the particle. 



1. Derivation of the mode- sum formula flS\ 
Consider the tail part of the scalar field, defined as 



$taii(aO = q Hm / G Iet [x, z(T')]dT' . (15) 



*0+ 



Recalling Eq. ([5]), we have 



F^{x) = qV a ^n(x) - q 2 (V a T^)G Ict [x,z(Tl)l (16) 



where G re t[a;, z(r_ )] = lim e ^ + G rc t[x, z(r_ — e)] contains only the smooth, tail part of the Green's function. Examine 
the second term on the RHS at the limit x — > z: By Eq. (33) of [5j, we have lim^z G Te t[x, z(tZ)] = R/12, where R 
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is the background scalar curvature. Since R — for Kerr, and since, as we mentioned already, the factor V q t_ is 
bounded at x — > z, we find that the second term on the RHS of Eq. (JTSJ) vanishes at this limit. Hence, 

F£a = gV Q $taii for x - z. (17) 
Next consider the direct part of the scalar field, defined as 



r-(x)+e 

$ dir = $ - $ tail = 9 Urn / G iet [x, z{T')]dT', (18) 

r_ (x) — e 



and define the direct /orce ^^(ie) replacing c — > (x)-l m ^1- ©• We nave 

i^z) = <7V Q $ dil .(x) + q 2 {V a T^)G, ct [x,z(TZ)], (19) 

since lim e ^ + G Ic t[x, z(r_ + e)] =0 by virtue of the retardation of the Green's function. Hence, for the direct field 
too, we find 

F^ = qV a <P dir fovx^z. (20) 

Finally, define the full force F a (x) by replacing J^ - ^ e — > J^-( x ' +e j n Eq, Q jpr, equivalently, through F a (x) = 
FtLi\( x ) + Fdir( x )]- Note F a (x) — qV a Q(x) holds precisely, for all x. Hence, defining 

F« s = gV Q $ rcs , F? = gV Q $ P , (21) 

we have, for all x, 

F t a ail (x) + F^(x) = F a (x) = F« B (x) + F?(x). (22) 

The goal of the somewhat elaborate construction in the last tew paragraphs is, partly, to establish the relations 
between the various 'forces' F(x) and their corresponding fields $(a;) — these will be needed in what follows. 
Based on Eq. (|2"2"|) we now write the scalar SF as 

F s a cU (z) = limits) 

X — >Z 

= lim[F x a es (x)-(F^(x)-F^(x))] 

X — >Z 

oo 

= lim Yu \f™W-(f™(*) -&"{*))]■ (23) 

m=— oo 

In the last step we have formally decomposed the tail force -F t aii( x ) ^--modes, introducing the notation 

f§T(x) = ± f F£(y, <p')e im ^-v'U<p', (24) 

where 'X' can stand for 'res', 'dir', or 'P', and y represents t,r,9. The m decomposition of F^ n (y,ip) is technically 
a standard Fourier series over the interval — n < ip < n (for any fixed y). Since F" ail {x) is continuous and piecewise 
differentiable over this interval, standard theorem of Fourier analysis (e.g., Theorem 4.2 of [35() guarantees that the 
sum over m in Eq. (f2U)l convergence to F^ lil (x), even at x = z. Moreover, the continuity of F t " ;i (x) assures that 
the Fourier sum converges uniformly (e.g., Theorem 4.4 of [HI), which allows us to switch the order of limit and 
summation in Eq. (I23| : 



OO 

(*) = E [It & - (/ d T (x) - ir o=))] • (25) 



m=0 



Here we have also folded the terms m < over onto m > 0, denoting = /? m + /v' m for to > 0, with 

fx' m ~° = /x' m_ °- Note that /x™ are real- valued, unlike f^ m which are complex (for m ^ 0). We stress that, ab 
initio, the convergence of the individual sums over / r "™, f£™, or /p m in Eq. (|23|) is not at all guaranteed. Also, one 
should not attempt to apply the limit in Eq. (|25|) to any of the three terms /™ individually. 
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We have now reached the crucial step of our derivation. In what follows we establish that, for any m > 0, 

C am = lim (/ d T(x) - = 0. (26) 

To show this, we start by inspecting the difference $dir — as x — > z. The explicit form of $dir was first worked out 
by Mino et al. [H, H3] ■ It can be written as [HI 

q qw(x) 

®dir(x) = —rr H 7 — \ — l~ const (for x near z), (27) 

e(x) e{x) 

where e = S 1 ^ 2 is the geodesic distance between x and the particle's worldline, and w is a certain regular function 
of x (and z), which has the form w — 0(Sx 2 ) near the particle. Recall Sx is short-hand for the coordinate difference 
5x a — x a — z a (t), where z a is a point on the worldline with the same t coordinate as x [hence, Sx a = (0, Sr, 59, Sip)]. 
To write the asymptotic form of $p near the particle, we first formally expand S near Sx = as S = So + S± + S2 + ■ ■ ■ , 
where S n represents the term of homogeneous order 0(Sx n+2 ). [The explicit form of Sq and Si was given in Eq. (J6j) 
above.] Starting from the definition of <&p in Eq. 1(7} we then have, for small Sx, 

- - (s^ky^ = e( i - km> + ° <Sxl) = ? + § + 0(fa2) - (28) 

From Eqs. ([UJ and igS} we obtain 

q (e 2 w - S2/2) /c 9N 

$dir-4'p = - L ^ 3 7 ; + const + 0{5x 2 ). (29) 

Note that both e 2 w and S2/2 are 0(<5x 4 ). Omitting the constant and the 0(Sx 2 ) term (neither can contribute to the 
eventual force at the limit Sx — > 0), we can therefore write the last expression in the form 

$dir-$P= "3 % (30) 

e o 

where eo = S*^ 2 , and Pmi is some multilinear function of the coordinate differences fa, of homogeneous order 0(fa 4 ). 
Recalling now Eqs. (|2H)) and (|2T|) . we finally obtain, at leading order in fa, 

n. - * - (i. - *p> - - f p " e ° v ° e ° , SW (31) 

where is yet another multilinear function of fa, of homogeneous order 0(fa 5 ). Note that $dir — ^p is continuous 
at the particle, and that Fg iT — Fp is bounded there, yet discontinuous (direction-dependent). 

Our goal, recall, is to confirm the validity of Eq. (|26|) . i.e., show that C am — 0. For this, we need to consider the 
in decomposition of Eq. (|3ip . at the limit x — * z. From Eqs. and ([2^)1 we obtain 
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where a m >o = 2, ao = 1, and <5y stands for 5r, 59 (recall St = 0). For later convenience, but without loss of generality, 
we have taken here the evaluation point z to be at ip = 0, so that Sip' = tp'. Note that, according to the definition of z, 
taking the limit x — > z results in taking x — > z, and hence also Sx — > (i.e., Sy — > as well as ip — > 0). Inspecting the 
integrand in the last equation, we remind that eg is a positive quadratic function of Sy and ip' which vanishes only at 
the particle, and Pj"] is a sum of a terms of the form a% n (p' k Sr n 59 5 ~ k ~ n , where fc, n are non-negative integers satisfying 
< k + n < 5, and a% n are constant coefficients (depending on z but not on r, 9). Below we analyze separately the 
terms k < 5 and k = 5, showing that both contributions to £ Qm vanish. 

Terms k < 5: We define R(Sy) = M~ 1 eo(S<p = 0), and, for fixed values of Sy (with small \5y\), split the integral 
in Eq. (|32|) into two domains: (i) \ip'\ < R and (ii) R < \ip'\ < ir. In domain (i) we can bound |P[f]| < c"P 5 and 
eo > C2P, where c" and C2 are some positive constants (depending on z but not on r,9). The absolute value of the 
integrand in Eq. (|32[) is thus bounded from above by c"/c|, and so the absolute value of the integral piece J_ R is 
< 2{c\/ c\)R. Since this vanishes at the limit fa — > (corresponding to i? — > 0), we conclude that the contribution 
to C am from domain (i) vanishes. Next consider domain (ii): Here we use eo > cal^'li |fa| < c aR- an d \S9\ < c$R 
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(where c 3 , a and c 5 are some other positive constants), to bound each of the above k,n terms of the integrand as 
\ a kn e o 5 f' k 5r n 69 5 - k - n cos[m(ip.-(f')]\ < c~ 5 cJc^ fe ~"|a£„</ £ - 5 i? 5 - fe |. The absolute value of the integral over domain 
(ii) is thus bounded from above by Cg \R 5 ~ k (n k ~ 4 — i? fe ~ 4 )| for k < 4 and by c" |i?log(7r/i?)| for k = 4, where Cg and 
c" are some positive constants. In both cases, the upper bound tends to zero as R —> 0, and we conclude that the 
domain (ii), too, gives a null contribution to C am . 

Term k = 5: The above upper-bound argument fails when k = 5. We now show that, nevertheless, the contribution 
to C am from this term vanishes as well, this time due to the symmetry of the integrand. For k — 5 (=>■ n = 0), the 
integral in Eq. (|32]) takes the form 

(¥>') 



^5,0 



■ cos[m(ip — ip')]dip'. (33) 



Once again, we split the integral into two domains, this time defined as (i) \<p'\ < R 1 / 2 and (ii) i? 1 / 2 < \ip'\ < tt. In 
domain (i), the absolute value of the integrand in Eq. (f33|) is bounded from above by , by virtue of e > C3 1 tp' | . 

Hence the absolute value of the integral piece is bounded from above by 2c~ 5 i? 1 / 2 , which vanishes as R — ► 0. 

Thus, the contribution to C am from domain (i) vanishes. Next consider domain (ii): For very small \5y\ and \ip\ we 
have here \ip'\ > R}l 2 > ISr/a] 1 ^ 2 ^> \6r\/c4, and similarly \ip'\ ^ \S8\/c 5 and \ip'\ 3> \f\. This allows us to expand 
the integrand in Eq. (|33j) about Sy, tp — 0, in the form 

' ~ ' ■ cos[m(^ - if/)] = cos(mtp') + 0(Sy, ip), (34) 



where we used eo(Sy = 0, tp') — Cg^ 2 |(^'|, eg being the coefficient of dip 2 in Cq. The integral of the leading-order term 
over domain (ii) vanishes from symmetry, since this term is anti-symmetric in tp' while the integration domain is 
symmetric. The integral over the 0(Sy, tp) terms clearly vanishes at the limit Sy, ip — » 0. We conclude, then, that the 
entire integral in Eq. (|33|) vanishes at the limit x —>■ z. 

We have shown that all contributions to C am vanish. Hence C am = 0. Importantly, the vanishing of C am does not 
depend on the direction from which the limit x — > z is taken. 

It is instructive to review the above analysis from the point of view of standard Fourier theory. The quantity C am 
is, formally, the Fourier transform of F^ h . — Fp (viewed at a function of <f for fixed y), summed over ±m for given m, 
and evaluated at x — > z. Eq. (|3ip implies that F^ h . — Fp is a pointwise continuous function of tp for any y (including 
y = 0). Standard Fourier theory [see, e.g., Theorem 4.2 of Ref. [35])] tells us that the Fourier series of such a function, 
evaluated at some point (p = ipo, converges to the average of the two one-sided values of the function at ipo (even if <pq 
is a step discontinuity). It is easy to show, based on Eq. ([3"Tj) . that the two-sided average of F^ ir — Fp (with respect 
to tp) vanishes at the limit Sx — > 0. This, reassuringly, is consistent with our finding C mn = 0. 

Having established Eq. ([26]) . we now return to our main line of development, and to Eq. (|25|) . We remind that, 
by construction, the limit x — > z in this equation is well defined and direction- independent. The vanishing of C mn 
therefore implies that each of the modes f" c T( x ) i s continuous: lim^^ fZT( x ) = fUTi 2 )- Our mode-sum expression 
for the tail part of the SF thus reduces to 



00 



m=0 

It now remains only to relate / r "™(z) to the numerical variables (/>™ s . From Eqs. (fTO)) and pTjl we obtain 
F t " s = q J2m=-oc ^/ a (4>7cs etmv )' where we have interchanged the derivative and the sum — this is allowed since 5> res 
is continuous and differentiable (see the discussion below) and hence its Fourier transform converges uniformly. On 
the other hand, we have F^ s = X)m=-oo /res™- Since both V Q (0™ s e lm ^) and / r "™ depend on <p only through a factor 
e mup j see (|24p], the orthogonality of the Fourier series implies 

f™ = < Z V Q (C s e imv ')- (36) 

In terms of the 'tilde' variables this becomes [recalling the definition of 0™ s in Eq. ([14]) ] f" c ™ =0 — (?V a </>™;r , and, for 
m > 0, 

= 2gV Q Re(e jmv> C s ) 

= ?V Q C S . (37) 

In the second equality we have made use of the symmetry relation <f>~™ — (</>™s)*i which is a direct consequence of 
$ rcs being a real field. Substituting f^{z) = gV Q 0™ s (z) in Eq. (I35|) finally establishes the mode-sum formula §\± 
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2. Convergence rate of the m-mode sum 

In practice, of course, one can only calculate a finite number of terms in the mode-sum series (|13[) . and so the 
question of convergence rate becomes important. To estimate the convergence rate of the series at large to, recall that 
the individual TO-mode contributions / r "" l (z) arise from taking the limit x — » z of continuous functions f^{x), which 
themselves are the Fourier components (multiplied by e lmv and summed over ±to for each to) of the function F^ as (x). 
The rate of convergence of the m-mode sum is therefore the rate of convergence of the Fourier series of F^ cs (x) at 
x — z. 

Here we touch upon a delicate issue: From Eq. (f2"2"|) we have 

F« s (x) = Ffa (x) + [F& (x) - F P « (x)} , (38) 

and we recall that F^(x) is continuous at z, but F^ I (x) — Fp(x) generally has a direction-dependence discontinuity 
there [see Eq. (|3"Tj) ]. Therefore, F£, s (x) is generally discontinuous at x = z, and, in general, it would have a jump 
discontinuity in the tp direction. Standard theorem in Fourier analysis (see, e.g., Theorem 4.2 of Ref. [351 ]) predicts, 
in this case, that the sum of Fourier modes converges very slowly: f^™(z) <~ 1/to at large to. Fortunately, however, 
the situation turns out better in our special case: As we established above, crucially, each m-mode of the difference 
Fjj- t — Fp happens to vanish at x = z (upon summation over ±m for each to) , which implies 

fr%?(z) = /»)• (39) 

In particular, we find that, at the point x = z, the convergence rate of the Fourier sum of F^ s (x) is the same as that 
of the Fourier sum of F^ a (x). The latter function is continuous and piece-wise differentiable at x — z (V^F t " u will 
generally have a direction-dependence discontinuity at that point). For such functions, we expect / t "™(z) ~ 1/to 2 at 
large to [see, e.g., Ref. [35J, Sec. 4.2.8]. Thus, we also expect /™™(z) ~ 1/m 2 , and expect that the mode-sum (fH{|) 
will, in general, converge like ~ 1/to. 



3. Behavior of 0™ s near the particle 

When devising a numerical integration scheme for the field equation (llip . it is important to know how regular 
the TO-mode fields 0™ s are near the particle. This, for example, influences the numerical convergence rate of a given 
finite-difference scheme. We therefore conclude our analysis by examining the regularity of </>™ s near the particle. 

From Eqs. (J5]) and (fT8j) we have $ res = <&taii + ($dir — 3>p). The tail field, as defined in Eq. (fl"5|) . is a continuous 
function of x near z. From Eq. (|29|) we have <£>dir — 3>p — > const as x — > z, and so this difference is continuous as well. 
Hence, $ res is continuous at x = z. It is then easy to show that the TO-modes 0™ s are also continuous at x = z. 

We next inspect the differentiability of 0™ s . From Eq. {57]), V Q 0™ S = f"™/q, where / r "™(x), as established above, 
are continuous at x = z. Hence 0™ s are differentiable at x = z. However, it is not the real fields </>™ s that interest us 
here, but rather the (complex-valued) fields 0™ s which would normally serve as variables for the numerical integration. 
To show that </>™ s , too, are differentiable at x = z, we first argue that the 'non-tilded' functions fj^ are continuous 
there, just like fj%P". To see this, consider the limit C am = linx,;^ (f£™ — fp m ), which can be evaluated in complete 
analogy to £ am , starting from Eq. (f3"2"]l with the replacements a m — > 1 and cos [to (ip—ip 1 )] — > exp[im(y>— tp')]. The same 
bounding arguments used for C am can apply for C am , merely replacing | cos[m((p—ip')]\ < 1 with | exp[im(<p— <p')]\ < 1. 
We find, once again, that only the term in Eq. (|34[) may potentially contribute to the integral at the limit x — > z; 
however, this time that contribution fails to vanish by symmetry — in fact, it yields 

nam i;„ i ram ram\ 5,0 / vP ) — imtfi' J / J 5/2 i OUQ TO, 

C = hm (/ dir -f P ) = — / e dp = < mW (40) 

x z L ~ J -it c 8 \(f | 5 [ 0, tor even to. 

[Note C a,m = —C a '~ m for all to, which, of course, is consistent with our previous finding, C am = 0.] Since the value 
of C am does not depend on the direction from which the limit x — * z is taken, we conclude that / di ™ — fp m is a 
continuous function of x near z. Recalling / r "™ = / t "™ + (f^ 1 — fp m ), and that / t "™ is continuous, we deduce that 
/ r "™ is continuous too. This, by virtue of Eq. (|36p . implies that the numerical m-mode variables 0™ s are differentiable 
at the particle. 

We note that the second derivatives of <^™ s are not necessarily continuous; indeed, they are not necessarily even 
bounded at the particle: We have gV /3 V Q $ rcs = V F T a cs = V p F^ n + V 3 (F d " ir - Fp), where V P F^ is bounded at 
the particle, but V' 3 (F^ ir — Fp) generally diverges there like ~ 1/eo [recall Eq. ((3T|) ]. Hence the second derivatives 
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V^V a $ res generally diverge at the particle, and their m-modes may, potentially, show a logarithmic divergence there 
(reminiscent of the logarithmic divergence of </> m (33|). 

In summary, when designing a numerical integration scheme for </>™ s , one should bare in mind that these fields are 
continuous and differ entiable (i.e., have continuous first derivatives), but not necessarily twice-differentiable. 



III. ELECTROMAGNETIC FIELD 



In this section we replace the scalar charge of the particle with an electric charge e, and formulate an m-mode-sum 
scheme for the EM SF acting on the particle. Once again, we assume that the particle is freely moving in a bound 
orbit around a Kerr black hole, which, neglecting SF effects, is a geodesic z(t) of the Kerr background. We will be 
interested in the value of the EM SF at an arbitrary point along this geodesic, and, as in the scalar case, we shall 
ignore the gravitational SF. The meaning of quantities like g a P , u", r, S, eo, and Sx a shall remain the same as in the 
scalar case, but 'Fx' and '/x" wu l now refer to the EM force (we avoid re-labeling the force as 'EM' in the interest of 
notational simplicity). The form of the m-mode-sum scheme in the EM case, and the analysis leading to it, resemble 
very closely those of the scalar case. This will allow us to skip many of the details of the derivation, pointing frequently 
to results from Sec. HU 



A. Preliminaries 



Let us denote the full (retarded) vector potential associated with our charge by A a , and assume that A a can be 
treated as a linear field over the fixed Kerr background. In the Lorenz gauge, the vector potential satisfies the sourced 
wave-like equation 

/oo 
S 4 [x ~ z{T)]{-g)- 1/2 u a {T)dT = -4irJ ai (41) 
-oo 

where J a is the charge-current density of the particle. Including SF effects to 0(e 2 ), the equation of motion of the 
particle reads pi H, Hf 

»uf> Vpu a = F s % = Urn Ffa (x) , (42) 

X — >Z 

where F s " lf is now the EM SF. In the EM case, the "tail force" field is given by 

F« a (x) = e 2 k a ^(x) lim f"^ * V 7 G$,[z, z(r')}u^' (r')dr' , (43) 

6 *® J — GO 

where the bi- vector G^, is the retarded Green's function associated with Eq. (141[) [with its un-primed (primed) indices 
corresponding to its first (second) argument] , V 7 acts on the first argument of G 1 ^, , and fc is a tensor defined at x 
through an extension of the velocity vector off the worldline: 

k a ^{x) = g a ^{x)u^{x) - g a(3 (x)u'(x). (44) 

The field u 7 (x) can be any smooth extension of it 7 ; the SF in Eq. (14"2")) depends only of the value of k on the worldline, 
and is not sensitive to the choice of extension. 

Eq. (|4"3"|) is analogous to Eq. (0J in Sec. [Til an d the same conclusions can be drawn in the EM case regarding the 
regularity of Ff^^x): This field is continuous and has at least piecewise continuous (bounded) derivatives. 



B. m-mode scheme: A prescription 

The procedure in the EM case follows very closely that of the scalar case. For a puncture function we now take[3| 

Al (s) - -4t faa + f IpUxSxt 3 ) , (45) 
ep(x) ' 

where ep is the same as in Eq. ([7]), and u a and are, respectively, the four velocity and the Kerr connection 
coefficients at point z(t). Recall, for given x, z(t) is the point along the particle's worldline intersected by the surface 
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t = const containing x, and Sx 13 = x^ — zP . One will note that, for small Sx, the parenthetical factor in Eq. (|45|) 
approximates the value of the four-velocity vector parallel propagated from z to x. We remind, however, that the 
definition in Eq. (|45p applies for any 5x, not necessarily small. 

As in the scalar case, we then define the 'residual' field A™ s , analogous to <& ros , through 

A™ = A a -Al, (46) 

and obtain a wave equation for A T ° S in the form 

V^V^ 03 = -4nJ a - V^VpAl = Z™. (47) 

We next formally decompose A™ s (t, r, 8, ip) into modes A r ° s ' rn (r, 8, ip)e tmv , just like in Eq. PH|) . and separate the ip 
dependence in Eq. (f47|) to obtain 

n (3)^ros,m = J_ f grcs e - t m V ' d(fi , _ Z rcs,m_ ( 48 ) 

Here Dy is a certain (m-dependent) second-order differential operator, which couples between the various vectorial 
components of A™ s > m , but not between different m-modes. Eq. ([48]) constitutes a set of four coupled hyperbolic 
equations, with a source term Z™ s ' m which, just like in the scalar case, extends away from the particle. Also in full 
analogy with the scalar case, the fields A r ^ s ' m are expected to be continuous and differentiable (have continuous first 
derivatives) on the worldline. 

Once again, we proceed by assuming that solutions A™ s > m to Eq. (|4"5|) , satisfying physical boundary conditions, 
have been obtained. The EM SF is then simply calculated through 

(49) 

x=z(t) 

where A r ^ s ' m (t, r, 8, ip) are real fields constructed from A r ^ s ' m (t, r, 8) through 

Ap m = 2Re (^ cs ' m e™ v ) for m > 0, and l^ es ' ro=0 = A™' m=0 . (50) 
As in the scalar case, the sum over m in Eq. (|4"9"|) is expected to converge at least as ~ 1/m. 



oo 

E 



C. Analysis 

The derivation of the mode-sum formula (149)) is entirely analogous to the derivation of Eq. (fl3)) in the scalar case. 
The tail part of the vector potential is given by 

Af xl (x) = e lim f ° G$,[x, z(r')]u ' (r')dT' , (51) 

which, recalling Eq. P^jL gives 

Ffaix) = ek a ^(x) [v,Af\x) - e(V 7 r_)G^[x,z(rr)]^'(rr)] . (52) 

Here G™p,[x, z(tZ)] = lim e ^o+ G z q^,[x,z(t- — e)] contains only the smooth, tail part of the Green's function, which, 
by Eq. (2.59) of Q, satisfies lim^^z GJj e J/[x, z(t_ )] = —Rppt/2 + gppiR/12. Since the background Ricci tensor Rpp 
and scalar curvature R both vanish in Kerr, and since, as mentioned above, the factor V q t_ is bounded at x — > z, 
we find that the second term in the square brackets in Eq. (|52[) vanishes at this limit. We hence obtain 

Ffa = ek^V^Af 1 forz^z, (53) 

which is analogous to Eq. (fT7|) . Following the same line of derivation as in Eqs. (p~8|) — (f2Tj) . we similarly obtain for the 
direct field F£ lt = ek^W^Af" (for x -> z), and for the full field F a = ek a ^\I 1 A p . The last equality holds for all x, 
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provided that the same off-worldline extension of k alSl is chosen for both F t " ;l and F^ ir — which we assume here. We 
then define, for all x, 

F T a es = ek a ^V~,Af s , Fg = ek a ^V 1 A p , (54) 

with fc"^ 7 extended the same way as in F a . 

With all quantities F£ defined as above for the EM case, we now reproduce all Eqs. ([22 |) -([25 |) . in their precise 
scalar-case form (with Fg, f^" 1 and f™ 1 now, of course, referring to the EM field). It should be noted that the 
value of each individual m-mode f^ m (or /™) will depend on the specific off-worldline extension chosen for A;"" 7 . 
However, the expression for the final SF in (the EM equivalent of) Eq. (|25|) will not depend on the fc-extension. 

As in the scalar case, the analysis proceeds by showing the validity of Eq. ([26]) . namely, that the difference f^ l — fp m 
vanishes at the limit x — ► z. We start with the expression for the asymptotic form of Ap hr (x) near z, as derived by 
Mino el al. [13] ■ We write it in the form [l2| 

, riir , . eupix) ewp(x) . 

Af T (x) = — P- 1 H P— - + const (for x near z). (55) 

c[x) e[x) 

where up is the four-velocity vector parallel propagated from z to x, and wp is a certain smooth function of x (and 
z) which has the local asymptotic form wp = 0(Sx 2 ), but whose exact form will not be important for us otherwise. 
For sufficiently small 5x we have up = up + Tp^uxSx 7 + 0(5x 2 ), where the remainder is a smooth function of x. 
Absorbing this remainder in the function wp, we rewrite Eq. (|55p as 

Af 1 = -(up+ Y^uxSx 1 ) + — ^ + const (for x near z). (56) 

Now, from the definition of the puncture function in Eq. (j45|) we obtain, recalling Eq. (|28|) . 



A? = - e (up+ f^fa 7 ) + + (Hfi.r> ,. C7) 

and so 



A dir_^p = ^ + const + Q(fa 2 ) 



e (e 2 wp - upS 2 /2) 
pf{5x) 



const + 0(Sx 2 ), (58) 

where Pi 4 ' is of homogeneous order 0{5x i ), and where we have absorbed the term arising from the difference between 
e and eo in the 0(5x 2 ) term. Using F^ iT — Fp = ek at31 V 7 (Ajf* — Ag) (for x — > z), we then obtain, at the limit x — > z, 

^r-^p= J V^» (59) 

where = ek af}l {z) (^effly-Pp^ — 3Pg 4 ^eoV 7 eoJ is 01 homogeneous order 0(5x 5 ). Eq. ([SI?)) in entirely analogous to 
Eq. (p{l"j) , and the analysis proceeds precisely along the lines of the discussion following the latter Equation, merely 
replacing PS, — > PS,. The conclusion, as in the scalar case, is that C am — 0, and, consequently, we deduce that Eq. 
(|55) is valid for the EM SF as well. 

We note that the validity of the crucial result C am = does not rely on any specific choice of the fc-extension, since 
only the value of k on the worldline enters Eq. (|59|) . Therefore, Eq. ([35)) is valid for any (smooth) extension of k. To 
proceed, we introduce a particular extension of k, denoted k, which we define through 

k a ^{x) ee k a ^{z). (60) 

Note fc a ^ 7 depends on t [through z{t)\ but not on r, 9, ip. We denote the functions P r " s and ffgp associated with the 
extension k by P r " s and / t "™ , respectively. Since the form of Eq. ([55)) is insensitive to the choice of extension, we may 
write 



oo 



m=0 
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where, just as in the scalar case, the 'tilde' notation indicates folding m < modes over to m > 0: / r "™ = /re^ + Zres m 
for m > 0, with j^ m =° = /£, m =°. 

It remains to relate the residual force modes f^™(z) to the numerical variables A r ^ s ' m . From Eq. (j54|) we have, 
choosing the k extension, F« s = e^_ M fc Q ^V 7 (^ es ' m e im *'). On the other hand, F r a cs = Em=-oo AT- In both 
sums, each of the m-terms depends on <p only through a factor e™ v , from which we deduce that 

/«™ = ek^V^A^e 1 ^). (62) 
In terms of the 'tilde' variables this becomes (recalling A r S s is a real field) 

j?™ = e P^V 7 4 cs ' ra (63) 

for all m > 0. Finally, taking the limit x — > z, we may remove the bar symbol from k, since all extensions yield the 
same value of k on the worldline. Substituting f" c ™(z) — efc a/37 (z)V 7 v4a es ' m (;z) in Eq. (|61[) establishes the mode-sum 
formula for the EM case, Eq. (149]) . 

The analysis of the convergence rate of the m-mode sum and the regularity of A T ^ s ' m near the particle is entirely 
analogous to the analysis in the scalar case, and we shall not reproduce it here. We state the main results: (i) The 
mode sum in Eq. (|49| is expected to converge at least like ~ 1/to. (ii) The variables A 1 p S,m are continuous and 
differentiable (i.e., have continuous first derivatives) on the worldline. (iii) The variables A r ^ s,m are not necessarily 
twice-differentiable on the particle's worldline. 

IV. GRAVITATIONAL PERTURBATIONS 

In this section we consider the gravitational SF acting on a test particle of mass /i, which is moving freely in a 
bound orbit around a Kerr black hole. At the limit fj, — ► the orbit is a geodesic of the Kerr background, again 
denoted z(t), and we shall prescribe an m-mode scheme analogous to the above for calculating the gravitational SF 
at an arbitrary point along this geodesic. Once again we retain the notation for 'Fx' and 'fx \ but these will now 
refer to the gravitational force. 

The physical nature of the gravitational SF is quite different from that of the EM (or scalar) forces: It is gauge 
dependent, and its physical interpretation is somewhat more subtle. This point is discussed in Ref. [38|, where, 
in particular, a general gauge-transformation law for the SF is derived. Our m-mode calculation scheme will be 
formulated within the Lorenz gauge (see below) — just like the fundamental formulation in Refs. Q and Q and like 
the standard mode-sum scheme. Practically, this means that the input for our new m-mode-sum formula will be the 
metric perturbation in the Lorenz gauge, and the output will be the Lorenz-gauge gravitational SF. The usefulness of 
the scheme stems from the recent advent of computational tools for direct integration of the perturbation equations 
in the Lorenz gauge [23l |39| . It should be stressed that the our mode-sum formula is not at all guaranteed to maintain 
its form in gauges other than Lorenz's. 

A. Preliminaries 

Let ha/3 denote the full (retarded) linear metric perturbation associated with our particle. Then define the "trace- 
reversed" perturbation, ^ a p = h a p — ^g a ph, where g a p is the background (Kerr) metric, and h = h a a is the trace 
of h a /3. We assume here that h a p is given in the Lorenz gauge, i.e., it satisfies the condition ^ a (i^ = 0. Then the 
perturbation satisfies 

/oc 
S A [x- z(T)}(-g)-^ 2 u a u p dT = -167rT Q/3 , (64) 
-co 

where T a p is the energy-momentum tensor associated with the particle, and R^a v p is the Riemann tensor of the Kerr 
background. The equation of motion of the particle, including 0(/i 2 ) gravitational SF effects, has precisely the same 
form as Eq. (|4*2"j) above, with the gravitational "tail force" field given by [1, 0] 

F t ln(x) = p?k a ^ s {x) lim r iX) \ s ^, 1 \x,z{T')}^\r')^'{T')dT'. (65) 

6 *^ J — oo 
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Here the bi-tensor GjS^, , is the retarded Green's function associated with Eq. (I64"]) . V5 acts on its first argument, 
and 

k a ^ s {x) = g a5 u p v?/2 - g al3 u"<u s - u a u^v?u 6 /2 + u a g^u s /A + 5 Q,5 / 7 /4, (66) 

where both g aS and m' 3 are evaluated at x, the latter through a smooth extension of the four-velocity off the worldline, 
just as in the EM case. In full analogy with the scalar and EM cases, F^ il (x) is continuous on the worldline, and has 
at least piecewise continuous (bounded) derivatives there. 

B. m-mode scheme: A prescription 

Introduce the puncture function 

= TJV\ + ( f ^ + f P-y Ua ) UxSx ^ > ( 67 ) 
where ep(ar), u a (t) and T„ 7 (i) are the same as in Eq. ([45]) . Then define the residual field through 

* r a e J = ~ (68) 

and obtain a wave equation for ^ x ^% in the form 

V^V 7 * r Q c | + 2i2" a " /J *~ - -lfaT afl - VW T <^ - 2^^ = Z-|. (69) 

Next formally decompose W£g(t,r,8,ip) into m-modes ^^(r, ( fi)e mip , just like in Eq. (|10p. and separate the <p 
dependence in Eq. l[6"9"|) to obtain 



n (3) , rcs,m _ J_ mea -imip' J J — ^res,m 

U G ~ 2-K J 13 ~ a/3 ' 



(70) 



where Dq is a certain second-order differential operator (depending on to), which couples between the various 
tensorial components of ip z ^g' m , but not between different m-modes. Equation (|7D|) constitutes a set of 10 coupled 
hyperbolic equations for the 10 components of V^'" 1 ' with an extended source term. The modes ip™p m are ' by the 
above construction, continuous and differentiable on the worldline. 

Assuming now that we have at hand solutions VVa'" 1 to Eq. (|70p . satisfying "physical" boundary conditions (in the 
sense discussed in Sec. [XXJ) , the gravitational SF is given by the simple formula 



F^[z(r)] = fi 



kaf 3jS J2 V5l/f S "' 



m=0 



(71) 



where if) a ' (t, r, 9, ip) are real fields constructed from the complex fields ipg ' (t, r, 6) through 



^ m = 2Re (f™> m e im *) for to > 0, and ^ m =° = ^ m=0 . (72) 

As in the scalar and EM cases, the sum over to in Eq. (|7Tj) is expected to converge at least as ~ 1/to. 

Below we explain the derivation of the mode-sum formula (|7"Tj) . referring to the scalar/EM cases for many of the 
details. 

C. Analysis 

We start with the formal expression for the tail part of the trace-reversed metric perturbation, 

*f»(x) = M lim f~" ' G$, 77 ,[x, z{r')y' {r')u^ {r')dT' . (73) 
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Comparing with Eq. (|65p we find 

Ffaix) = »k^ s (x) [v,*$ l (z) - »(V S T-)G$,^,[x, z{rZ)V ' (rZ)u< ' (rZ)] , (74) 

where G^p, ,[x, z(tZ)] = lim e ^ + G 1 ^, , [x, z(j- — e)] contains only the smooth, tail part of the Green's function. 
Unlike in the scalar and EM cases, here the factor oc G ret generally does not vanish at the limit x — ► z: From Eq. 
(2.11) of [3| we have lim x _^ z G^ C J, 77 , [x, z(tZ)] = Rp^piy + Rpyp'-y. However, it may be readily verified that the 

contraction k af3lS Gp c p,^,u t3 u 1 in Eq. ([71)1 does vanish at x — > z, for any value of a and S. (Some of the terms in 
this contraction vanish by virtue of the symmetry of the Ricmann tensor, and the others due to the Ricci-flatness of 
the Kerr spacetime). Thus, 

F« a = /xfc^V^ 1 toix^z, (75) 

in analogy with Eqs. fT^ and ([53|). We similarly obtain F d Q h . = fik af3 ~« 5 \7 s i$f* (for x -> z), and F a = fik^^Vs^p-y 
(for all a;), and for all x define F T a cs = fik a ^ s \7 s ^ and Fp = fik al3 ~< s V 's^%, with the same off-worldline extension 

of k a ^ s chosen in all cases. Once again we write the SF as in Eq. (f2"5"]) . and proceed to show that the limit C am , 
defined in Eq. (|26)l . is null in the gravitational case too. 

The form of the direct part of the Lorenz-gauge perturbation was derived by Mino el al. [13] ■ It can be written as 


. = + + congt or x near 

' ' e(x) e(x) 

where wp 7 is a certain smooth function of x (and z) which has the local asymptotic form — 0(Sx 2 ). [tip, recall, 
is the four- velocity parallel propagated from z(t) to x\. As in the EM case, we replace tip — > up + Tp S ti\5x s in Eq. 
([76]) . absorbing the (smooth) 0(8x 2 ) error within the function wp 7 . We obtain 

Vf* = — [upu 1 + (fp 5 u 7 + T^gup) u\5x s ] + M Wpi + const (for x near z). (77) 



We wish to compare this local asymptotic form with that of the puncture function defined in Eq. (|67|) . The latter 
reads, in analogy with Eq. ([57|) of the EM case, 



= ^ [upti, + (f ^ 7 + f^up) u x Sx s ] + g^gfcgg + (Sx 2 ). (78) 

Thus, *^-*^ 7 = [ie~ 3 (e 2 wp 7 - upu- f S 2 /2) +const+0(fe 2 ), which has the form eZ 3 Pp^ (5 x)+ const + 0(6 x 2 ), where 

is a certain smooth function of a; (and z), of homogeneous order 0(5x i ). Using Fjj- T — Fp = [ik a P lS V s{'^'pl l — ^Jj 7 ) 
(for x — ► z), we obtain, at the limit x — > z, 

P<% (6x) 

F^~Ft= J ^, (79) 
e o 

where P^ = ek af3j5 (z) (^VsPp^ — SPp^toVgeoj is a smooth function of x (and z), of homogeneous order 0(5x 5 ). 
This, once again, has precisely the form of the equivalent scalar-case expression, Eq. ([31]) , and a similar analysis (merely 
replacing P^ — > P&) then leads to the conclusion C am — 0. We find that Eq. (|35|) is valid in the gravitational case 
too. 

The derivation proceeds just as in the EM case: We introduce the "fixed contravariant components" /c-extension, 
k a ^ s (x) = k a ^ s (z), for which f™ = ^ik a ^ s V s (^' m e tmv ). In terms of the 'tilde' variables defined in Eq. $72]) 

this becomes (recalling 4^ is a real field) = fik al3 "> s V s ^' m and so f?™(z) = fik a/3 ^ ; (z)V 'sftp*'™ '(z) , where 
we have removed the bar off k since all extensions coincide at x = z. Finally, choosing the k extension in Eq. (|35|) 
and substituting for f^™(z) from the last expression, we arrive at the mode-sum formula (|7Tj) . 

Once again, the analysis of the convergence rate of the m-mode sum and the regularity of the residual field replicates 
the scalar-case analysis of Sec. UIl so we merely re-state the results for the gravitational case: (i) The mode sum in 
Eq. (|7ip is expected to converge at least like ~ 1/m. (ii) The variables ipgL' m are continuous and differ entiable on 
the worldline, but are not necessarily twice-differentiable there. 
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V. SUMMARY AND CONCLUDING REMARKS 



Equations (fT3|) . (|49|) and ([7T|) prescribe the construction of the scalar, EM, and gravitational SFs, respectively, 
within our new m-mode-sum scheme. In each case, the raw input for the SF formula are the m-modes of the residual 
function (4>™ s , A r ^ s,m , or ip^ ,m ), which are to be obtained through numerical time-evolution of a 'punctured' version 
of the field equations. Once these residual modes are calculated, they require no further regularization: The SF is 
given as a simple sum of (certain combinations of) the derivatives of these modes, evaluated at the particle's location. 

The most computationally-involved stage in the implementation of the scheme is, of course, the calculation of the 
?n-modcs of the residual function, which requires numerical evolution in 2+1D. The feasibility of such calculations was 
demonstrated in Ref. [33| for the test case of a scalar field in Schwarzschild. Since the numerical method developed 
in Ref. [331 ] does not rely on the spherical symmetry of the background spacetime, it is directly applicable to a Kerr 
background (work to extend the existing scalar code to the Kerr case is underway [1(3]). We envisage applying a 
similar numerical method to evolve the coupled set of Lorenz-gauge metric perturbations, but this will require further 
development. 

Working in 2+1D, and calculating the SF using the new m-mode scheme, would offer significant practical advan- 
tages, especially in the Kerr case. First, the standard ^-mode regularization procedure involves the projection of the 
spheroidal-harmonics components into a basis of spherical harmonics, which greatly complicates the calculation. This 
complication is spared within the new m-mode scheme, which avoids the multipole decomposition altogether. Second, 
in the standard £-mode scheme the numerical variables (i.e, the £, m modes of the full fields) require regularization, 
which reduces their numerical accuracy: the full modes must be calculated with great initial accuracy in order for the 
final, regularized SF to be only moderately accurate (this point is elaborated on in, e.g., [23|). In contrast, the numer- 
ical variables in the m-mode scheme (i.e, the m-modes of the residual fields) require no regularization, which loosens 
the numerical accuracy requirements. Third, in the i'-mode scheme one has to analyze separately the contribution to 
the mode sum from the "non-radiative" low multipoles (t = 0, 1 in the Schwarzschild case); in the Kerr case, merely 
identifying this contribution is not a trivial task. The new m-mode scheme, on the other hand, requires no special 
treatment of the low multipoles — their contribution is automatically contained within each of the m-modes. 

This work focuses on establishing the theoretical grounds for the m-mode scheme, and it does not discuss the 
practicality of the numerical implementation. Many of these practical issues are discussed in Ref. 33]. Here we just 
mention one such implementation issue, which concerns the freedom in choosing the puncture function. The particular 
form selected for <E>p in Eq. ([7]) [and for and in Eqs. (jl5]) and (f6"T|) . respectively] is, of course, not unique. 
Focusing, for the moment, on the scalar case, consider the class of puncture functions 

$ P = q =, (80) 

VS + Si + AS' v ; 

where AS is any smooth function of the coordinates which vanishes at the particle at least as 0(Sx 4 ). Now imagine 
that we re- formulate the m-mode scheme by replacing $p — * <I>p (for a specific, but arbitrary, choice of AS). The 
analysis in Sec. Ill CI can then be repeated step by step, with only a slight adjustment: In Eqs. (f28|) and ([29]) . we 
should replace S2 — > S2 ~ AS. Since S2 — AS, just like S%, is a smooth function of order 0(Sx 4 ), Eq. (|2"5j) will 
maintain its form, with only the explicit form of Pmi (Sx) being affected. The rest of the analysis does not depend 
in any way on the explicit form of Pw . In particular, we find that the final mode-sum formula (| 13() is applicable for 
any choice of a puncture function within the class <I>p. In just the same way, we can generalize the puncture schemes 
in the EM and gravitational cases too, by introducing generalized classes of puncture functions obtained by taking 
ep(x) = V<So + Si + AS in Eqs. (|45|) and (|67|) . So long as AS is a smooth function of the coordinates which vanishes 
at the particle at least as 0(Sx 4 ), the form of the mode-sum formulas in Eqs. and (|TTj) will not change. 

The above freedom in choosing AS* can be exploited in the actual implementation of the puncture scheme, in order to 
simplify the Fourier integrations involved. In Ref. [33j, for instance, we have made the replacement Sip 2 — > 2(1 — cos Sip) 
in the expression for So (which amounts, in the case considered there, to specifying a certain non-zero function AS*), 
and this allowed us to calculate analytically all the necessary Fourier integrals. Other choices may simplify the 
implementation in other cases. One should be careful, though, to avoid choices of AS which nullify So + Si + AS at 
points other than the particle's location. 

Finally, we comment that the theoretical foundation for our new m-mode scheme, herein established based on 
the 'tail' + 'direct' decomposition, can alternatively be formulated based on Detweiler and Whiting's 'R' + 'S' 
decomposition [6j. The mathematical details of the analysis would be quite similar, leading, of course, to the same 
final m-mode formula for the SF. 
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Our residual field <E> rcs is not to be confused with the "R- field" introduced by Detweiler and Whiting [6J: The latter, by 
definition, is a solution of the vacuum field equation, whereas our function $ rcs is not necessarily so. 

In this regard, we point to a mistake in Eq. (1) of Re f. II ill : This form of the scalar-particle equation of motion is, in 
general, not self-consistent, since the SF, as defined in [llj, is not necessarily co-alighned with the four-acceleration. The 
correct form is given in Eq. @ here. 

Throughout this work we take the principal values of the Boyer-Lindquist azimuthal coordinate to lie in the range — n < 

if < 7T. 

In what follows we use the label 'P' alternately as upper- or lower-script, and similarly with the labels 'res', 'tail', and 
'dir'. This should not cause any confusion. 



